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Abstract 

We investigate a class of power series occurring in some problems in 
quantum optics. Their coefficients are either Gegenbauer or Laguerre 
polynomials multiplied by binomial coefficients. Although their sums 
have been known for a long time, we employ here a different method 
to recover them as higher-order derivatives of the generating function of 
the given orthogonal polynomials. The key point in our proof consists in 
exploiting a specific functional equation satisfied by the generating func- 
tion in conjunction with Cauchy's integral formula for the derivatives of 
a holomorphic function. Special or limiting cases of Gegenbauer poly- 
nomials include the Legendre and Chebyshev polynomials. The series of 
Hermite polynomials is treated in a straightforward way, as well as an 
asymptotic case of either the Gegenbauer or the Laguerre series. Fur- 
ther, we have succeeded in evaluating the sum of a similar power series 
which is a higher-order derivative of Mehler's generating function. As a 
prerequisite, we have used a convenient factorization of the latter that 
enabled us to employ a particular Laguerre expansion. Mehler's summa- 
tion formula is then applied in quantum mechanics in order to retrieve 
the propagator of a linear harmonic oscillator. 

1 INTRODUCTION 

The problem of evaluating the Zth-order correlation function ((a') 1 a 1 ) for 



1 



a one-mode squeezed thermal state of the quantum electromagnetic field by 
employing the photon-number distribution led one of us [Tj to consider the 
following power series involving Legendre polynomials: 



G] 2) (t,w) := ( ? ) t n P n (w), (w:=cos9, \t\ < 1, I = 0, 1, 2, 3, ...). 




introduce the square root of a complex variable as follows: 



q : = (1 - 2wi + t 2 )5, 



(i = 




(1.2) 



The sum of any series of ascending powers (jl-l|) is the product of three factors: 
the generating function of the Legendre polynomials, the Zth power of a specific 
variable, and the Legendre polynomial of degree / of another specific variable: 



Note that the series expansion (jl.ip reduces for Z = to the power series of 
the generating function of the Legendre polynomials. 

A similar result has been found later for the Laguerre polynomials [2]. 
We obtained a formula which is useful when the general solution of the quan- 
tum optical master equation describing dissipation is applied to a single-mode 
displaced squeezed thermal state [3]- 

When we performed the above-mentioned summations we were not aware 
of the fact that the corresponding formulae and the similar ones for the Gegen- 
bauer and Hermite polynomials have been known for a long time [3]. This 
happened because such formulae are not mentioned in any book devoted to 
special functions except for the cited textbook. Rainville derived them by 
making combined use of an appropriate double power series and of the gen- 
erating function of the orthogonal polynomials under discussion. Owing to 
their importance for some physical problems, we find it useful to give in this 
work a comprehensive overview based on an alternative approach. We prove 
that a property of the type (|1.3p is shared by two one-parameter families of 
classical orthogonal polynomials, namely, the Gegenbauer and Laguerre poly- 
nomials. The main ingredients of our method are Cauchy's integral formula 
for the derivatives of a holomorphic function and specific functional equations 
satisfied by the generating functions of these well-known polynomials. 

The outline of the article is as follows. Section 2 deals with the Gegen- 
bauer series, while in Sec. 3 the obtained formula is applied to the particular 
cases of the Legendre and Chebyshev polynomials. In Sec. 4 we review the 





(1.3) 
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Laguerre series. The sum of the series of Hermite polynomials is recovered in 
Sec. 5 in a straightforward manner, as well as an asymptotic limit of both 
the Gegenbauer and Laguerre series. In Sec. 6 we compute a similar power 
series of an arbitrary derivative of Mehler's generating function, finding a new 
formula. By way of illustration, Mehler's summation formula is employed in 
Sec. 7 to write in closed form the well-known propagator of a linear harmonic 
oscillator. Concluding remarks are presented in Sec. 8. 

2 GEGENBAUER POLYNOMIALS 

The generating function of the Gegenbauer polynomials of order A, 
G { X) (t,w) := (1 -2wt + t 2 Y x , (Ay^O, A > --, w:=cos0), (2.1) 
is holomorphic in the variable t for \t\ < 1 and has the Taylor expansion [5 J 

oo 

G [X) (t,w) = J2 tn °n(w), (w := cos 9, \t\ < 1). (2.2) 

n=0 

Accordingly, Cauchy's integral formula for the derivatives of a holomorphic 
function (CIFD) [BJ[7] gives the following integral representation of the Gegen- 
bauer polynomial C x (w) of degree I and order A: 

> 1 f {0+) 1 (X) 

C * {w) = 2riJ * (2.3) 

The contour of integration in eq. (|2.3p is a simple loop encircling the origin 
t = counterclockwise and situated entirely within the open disk \t\ < 1. Note 
that the sum of the power series of interest, 

G{ A) (i, W ):=JV " )t n C*(w), (w=:cos9, \t\ < 1, I = 0, 1, 2, 3, ...), 

(2-4) 

involves the £th-order derivative of the generating function: 

G\ X \t,w) = ^^ l Gi X \t,w). (2.5) 

We employ CIFD again and find, after a translation of the variable of integra- 
tion, 

G\ X \t,w) = ^- j 0+) ds-^G { X) (s + t,w). (2.6) 
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We now take advantage of the following functional equation for the generating 
function ([23]) : 

G { X) (s + t,w) = G {X) (t, w) G ( X) (t, r,), (2.7) 



where we have denoted 



s w — t 
£:=-, 77:= . (2i 



Insertion of eq. (|2.7p into eq. (j2.6|) . followed by the change of the integration 
variable to £ defined in eq. (|2.8[) . yields 



77). (2.9) 



G w (t, «,) = g?) ( t, «,) (i) ' JL / (0+) * 4 A) (e, 

For real t, the obvious inequality r/ 2 ^ 1 allows us to use eqs. (|2.ip and ()2.3p 
and find the formula 

(A/0, ^>-^> w:=cos0, |t|<l, i = 0,1,2,3,...). (2.10) 

The validity of eq. (|2,10p can be extended by analytic continuation to any 
complex value t inside the unit circle. An equivalent form of eq. ()2.10p . 

t('T) rcL.™ = <r 2A ^ (n^i) , (2.ii) 

m=0 ^ ' \ 9 / 

coincides with the formula written by Rainville [8 : . 

3 LEGENDRE AND CHEBYSHEV 
POLYNOMIALS 



Recall that the Legendre polynomials are the Gegenbauer polynomials 
of order A = \: 

P l {w) = Cf{w). (3.1) 

Therefore, by setting A = \ in eqs. (|2.4p and (|2.10p . we recover eqs. (jl.ip and 
(|1.3p . respectively. For 
chosen by Rainville [9]: 

±(' + '")fP l+m{w ) = ^) 1+ 'p l (^.). (3.2) 
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(jl.3p . respectively. For A = ^, eq. (|2.1ip reduces to the equivalent formula 



m=0 
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The Chebyshev polynomials of the first kind, 

Ti(w) := cos(W), {w:=cos6, 1 = 0,1,2,3,...), (3.3) 
are the limiting case A = of the Gegenbauer polynomials, as follows [10] : 

ToH = C », (3.4) 



^H = ^Q°H, (/ = 1,2,3,...), 



with [IT] 



C » := 1, 



C?W := lim 

A— >0 



(1 = 1,2,3,...). 



Hence the generating function of the polynomials (j3.7|) is the derivative 



G, 



(0) 



(3.5) 

(3.6) 
(3.7) 

(3.8) 



A=0 



while the functions studied here are given by the limit 



G^iw) := lim 

' A^O 



1 r W 



Gr(t,w) 



(1 = 1,2,3,...). 



(3.9) 



The generating function (j3.8f) . whose explicit form is the principal de- 
termination of the logarithm 



G { o\t,w) = -ln(l-2wt + t 2 ), 



is the sum of the power series [12 



G$\t,w) = J2t n± T n (w) 



(w := cos 8, \t\ < 1). 



(3.10) 



(3.11) 



n=l 



In turn, eq. (|2.4|) leads us to write a function (|3.9|) as the sum of the series 
G| 0) (t,^) = £( " )t n -T n (w), (w:=cos6, \t\ < 1, / = 1, 2, 3, ...), 

71=i 

(3.12) 



while eq. f)2.10[) gives its explicit expression: 



G^^w) 



t\ 2 /'lo-f 



qj I 



{1 = 1,2,3, 



(3.13) 



In particular, for I = 1 |13j . 

2wt + 1 2 ' 



Y t n T n (w) = 1 Wt (3.14) 
^ K 1 l-2wt + t 2 

n=0 



Equations (|3.12p and ()3.13p also read: 



(3.15) 



The Chebyshev polynomials of the second kind, 



Ui(w) := Sm[ ^0 )d \ (w:=cose, I = 0, 1, 2, 3, ...), (3.16) 



are the Gegenbauer polynomials of order A = 1 

U l {w) = Cl{w). (3.17) 
Accordingly, the sum of the series 



Gf\t,w) :=JV U At n U n (w), 

n=l ^ ' 



is the special case A = 1 of eq. (|2.10p : 



(w:=cos0, \t\ < 1, 1 = 0,1,2,3,...) 

(3.18) 



Gf\t,w) = q- 2 (^jUi (^-^\ , (1 = 0,1,2,3,...). (3.19) 

An equivalent formula is given by eq. (|2.11|) written for A = 1: 
i + m \ , mr . , . / 1 \ TT w — t 



m=0 



4 LAGUERRE POLYNOMIALS 

For the sake of completeness, we prove the similar property of the La- 
merre polynomials, repeating step by step the argument from Section 2. The 
generating function of the Laguerre polynomials of order a, 



(a > -1), (4.1) 
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is holomorphic in the variable t inside the circle \t\ = 1 and has the Taylor 
expansion |15j 

oo 

si a \t,u) = Y,t n K(u), (1*1 <1). (4.2) 

n=0 

Making use of CIFD, we write the Laguerre polynomial Lf(u) of degree I and 
order a as an integral representation [16] : 

L?{U) = 2^J dt Wl S { a \t,u). (4.3) 

Just as in eq. (I2.3p . the contour of integration encircles the origin t = 
once, counterclockwise, inside the unit circle. Now the power series under 
investigation, 



(u), (|t|<l, Z = 0,1,2,3,...), (4.4) 



n=l 

has the sum 



5?° ) (*,«) = |(^) , 5W(t,tt). (4.5) 

Applying once again CIFD, we get with a translation of the variable of inte- 
gration: 

S > ' (t ' u) = 2^7 d S7ri Si a \ S + t,u). (4.6) 
The generating function (|4,ip satisfies the functional equation 

Combining it with the change of variable 



(4.7) 



in the integral, eq. (|4.6p becomes 

S <«»( t ,„) = #><«,„) (_!_)' _L/ <0+, rf(? L5r (C, ^) • (4.9, 

On account of eqs. ()4.ip and (|4.3p . eq. (|4.9p finally yields the sum of the series 
(g3D: 

(a>-l, |t| < 1, 1 = 0,1,2,3,...). (4.10) 
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Note that eqs. (|4.4p and (|4.10p are equivalent to the formula preferred by 
Rainville [TTJ: 

t ( 1 T ) ^ (-r- exp 1 If (^) . 

(4.11) 

5 HERMITE POLYNOMIALS 



The generating function of the Hermite polynomials is the exponential 
W (z, x) := exp(2xz - z 2 ). (5.1) 
Being an entire function, the Taylor series |18] 



W Q (z,x) = y2 — H n (x) (5.2) 



n=0 



converges in the whole complex z-plane. The more general series we are in- 
terested in, 

Wi(z,x) :=J3( " \— y H n {x), (Z = 0,1,2,3,...), (5.3) 

n=l ^ ' n ' 

has the sum 

W l (z,x)) = ^(jA W (z,x). (5.4) 

In order to evaluate it, we apply the method of Sees. 2 and 3, taking advantage 
of the functional equation satisfied by the generating function (|5.ip : 



W (v + z, x) = W (z, x) W (v, x - z). (5.5) 

The result is 

t l 

W l (z,x) = exp(2xz-z 2 )-H l (x-z), (1 = 0,1,2,3,...), (5.6) 
and may be written in the equivalent form chosen by Rainville |19j : 

V" — Hi +m (x) = ex.p(2xz - z 2 ) Hi(x - z). (5.7) 
^— ' ml 

m=0 

However, eq. (|5.6|) can be established as an asymptotic case of the similar 
formula either for the Gegenbauer polynomials, eq. (|2.10p . or for the Laguerre 
polynomials, eq. (j4. 10p . 
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In the first case, our starting point is the limit |20| 



Hence 



W (z,x)=1im G ( X) {\-h,\-^x) 

A— >oo 



\:Hi{x) = lim \\-^Cf{\-^x) 

l\ A— s-oo L 



Making use of eqs. (|5.4|) and (j2.5|) . we get 



Wi(z,x)= lim G^(A ?z,\ 2x) 

A— >oo 



(5.8) 
(5.9) 

(5.10) 



Substitution of eqs. (|2.10p . (|5.8p . and ()5.9[) into eq. (|5.10p yields again the 
result flUSH- 

In the second case, we start from the identity [21] 



W (z,x) = lim { S$ 



(a) 



2 \ 2 

— I z, a 

a 



which gives the limit 
1 



U,(.r) lim I I -V IS* ( a 



2\ 5 



o 



(5.11) 



(5.12) 



On the other hand, taking note of eqs. (|5.4p and (|4.5p . we find: 



Wj(«,x) = lim ^ 5 



(a) / / 2 \ 2 



— I z, a 

a 



i - I -)x 

a 



(5.13) 



By inserting eq. (|4.1(J|) as well as the limits (|5.11|) and (|5.12|) into eq. (|5.13|) . 
we recover once more eq. (j5.6|) . 



6 A GENERALIZATION OF MEHLER'S 
SUMMATION FORMULA 



We intend to evaluate the sum of the power series 

OO / v n 

M^z^y) :=J2( ? J ZJ7^. H n (x)H n (y), (\z\<l. 1 = 0. J . 2. .. 



I J n!2 

Its special case £ = is Mehler's power series expansion 



z 

M (z, x,y) = Y, ^ H n {x)H n (y), (\z\ < 1), 



n=0 



(6.1) 
(6.2) 
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of the generating function \22\ \ 

M {z,x,y) := (l - z 2 p exp j _L_ [ 2 xyz - (x 2 + y 2 ) z 2 ] \ . (6.3) 

The analytic function ()6.3p factors in the following way |25| : 

M Q (z, x, y) = S H) (z, \{x- y) 2 ) S H) (-z, \{x + y) 2 ) , (\z\ < 1), 

(6.4) 

where Sq 2 (t, u) is the generating function (|4.ip of the Laguerre polynomials 
of order a = — |. On the other hand, the function f)6. 1 1) involves a /th-order 
partial derivative of Mehler's generating function (j6.3|) : 



/ 5 \ 1 

M l (z,x,y)) = -(—J M (z,x,y). (6.5) 

Insertion of the factorization (|6.4p into eq. (|6.5p . followed by application of 
Leibniz' rule and use of eq. (|4.5p with a = — 5, leads us to a finite expansion: 



M l (z,x,y)= [z,l(x-y) 2 )st h) (-z^x + y) 2 

{\z\<l, I = 0,1,2,3,...). (6.6) 

By substituting the series sum (|4.10p and employing again eq. ()6.4p . we estab- 
lish the following formula: 

Mi(z,x,y) = M o(z,x,y) x 



m=0 



1 + zJ l - m \2{l-z)J \2(l + z 



(\z\<l, 1 = 0,1, 2,3,...). (6.7) 

The function ()6.7p is a finite sum involving Laguerre polynomials of order 
a = — 7^. It can equally be expressed in terms of Hermite polynomials, via the 
identity [26] 

H 2m (x) = {-\) m m\ 2 2m LZ?(x 2 ). (6.8) 
We get thus an alternative explicit formula: 



-1)' / z 

Mi(z,x,y) = [yZT^ ) M o( z , x ,v) x 

1 



m=0 



I \ f 1 — z\ ( x — y \ I x ~\- y 
m ) {-—z ) H ^- m) I ~ 1 H2m ' — 



[2(1 -z)] 2; V[ 2 (1+^X ; 



(|z|<l, / = 0,1,2,3,...). (6.9) 
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7 PROPAGATOR OF A LINEAR HARMONIC 

OSCILLATOR 



In nonrelativistic quantum mechanics it is instructive to evaluate the 
propagator of a linear harmonic oscillator by employing Mehler's summation 
formula, eqs. (|6.2|) - (|6.3p . In the coordinate representation, the propagator is 
defined as the probability amplitude for finding the particle at some point x b , 
at time t b , when it has originally been at another point x a , at time t a : 



K(x b , t bl x a , t a ) . — {x b: t b | x a , t a ) . 



(7.1) 



For an oscillator of mass m and classical angular frequency ui it is convenient 
to denote a := (^) 5 > 0. We recall the eigenvalues 



n + -]Twj, 



(n = 0,1,2,3,...), 



(7.2) 



as well as the corresponding eigenfunctions of its energy operator H, 



u n (x) 



a 



ir l l 2 n\2 r 



1/2 



cxp 



1 



a 2 x 2 ) H n (ax), (n = 0, 1,2, 3, 



Making use of the time-evolution operator in the Schrodinger picture, 



-r(tb-t a )H 
n 



U{t b , t a ) = exp 

the transition probability amplitude (|7.ip reads: 

(x b ,t b | x a ,t a ) = (x b \U(tb, t a )\x a ). 



(7.3) 
(7.4) 

(7.5) 



Insertion of eq. (|7.4p into eq. (|7.5p yields, via the spectral decomposition of 
the energy H, the standard eigenfunction expansion of the propagator: 



K(x b ,t b ;x a ,t a ) = y^exp 



71=0 



E 
n 



u n (x b )u* n (x a ). (7.6) 



Substitution of the energy levels (|7.2jl and eigenfunctions (|7.3p into eq. ([7.6 
enables us to exploit Mehler's expansion (|6.2|) to get a compact formula: 



K 



. . a 

[X b ,t b ] x aj t a ) = —< 



exp 



1 



-i^{t b -t a ) - - a (x b + x a ) 



M Q (exp [-iuj(t b - t a )] , ax b , ax a ) . 



(7.7) 
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Equations (|6.3p and (|7.7p finally give the exact closed-form expression of the 
propagator of the one-dimensional harmonic oscillator [28] : 

K(Xb,tb',X a ,t a ) = < — — ; : T- — tt > X 

l_27rftsin[w(t 6 - t a )\ J 

ex P f oh ■ iwT — Tvi { ( x b + x a) cos [w(t 6 - t a )] - 2x b Xa] ) . (7.8) 
\2hsm[uj{t b - t a )\ 'J 

Since Mehler's summation of the eigenfunction series (|7.6[) does not always 
fulfil the convergence requirement from eq. (|6.2|) . the propagator (|7.8|) is a 
distribution involving 5-type singularities when the time difference t b — t a is a 
multiple of the classical half-period of oscillation tt/uj. For instance, eq. (|7.5p . 
and eq. (|7.6p as well, displays the equal-time values 

K(x b ,t;x a ,t) = S(x b - x a ). (7.9) 

We stress that the transition probability amplitudes of the type (|7.ip 
are at heart of Feynman's path-integral formulation of quantum mechanics. 
In particular, for the linear harmonic oscillator the propagator (|7.ip has been 
evaluated exactly in this framework |29[ 130] . However, the explicit formula 
(|7.8p can also be obtained in an alternative simpler way [28]. Had one started 
from eq. (|7.8|) . one could then find, conversely, the energy levels (|7.2p and 
eigenfunctions (|7.3p . Indeed, insertion of Mehler's power series (|6,2p into 
eq. (|7.7p allows one to compare the resulting formula with the general eigen- 
function expansion (|7.6p . One thus readily gets the well-known solution of 
this energy eigenvalue problem. 

8 CONCLUSIONS 

The sum of a series of the type (|1.1|) involves the /th-order derivative of 
the generating function of the orthogonal polynomials. We have evaluated it 
in a straightforward way for both the Gegenbauer and Laguerre polynomials 
taking advantage of the analytic structure of their generating functions. The 
explicit sum is found to be a product of three factors: the generating function, 
the Zth power of a specific variable, and the corresponding orthogonal polyno- 
mial of degree I in another specific variable. However, the limiting case of the 
Chebyshev polynomials of the first kind stands out since their generating func- 
tion is missing in the explicit formula. Note that the above-mentioned result 
is by no means general: it holds neither for arbitrary Jacobi polynomials, nor 
for some nonclassical orthogonal polynomials [27]. This is due to the analytic 
peculiarities of the generating functions in question. Along the same lines, we 
have succeeded in evaluating the higher-order derivatives of Mehler's generat- 
ing function as finite sums involving products of pairs of Hermite polynomials 
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of even degrees, eq. (|6.9p . The importance of Mehler's formula in quantum 
mechanics is finally emphasized. 
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